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Glasses have a large excess of low-frequency vibrational modes in comparison with continuous 
elastic body, the so-called Boson Peak, which appears to correlate with several crucial properties of 
glasses, such as transport or fragility. I review recent resultsii^ showing that the Boson Peak is a 
necessary consequence of the weak connectivity of the solid. I explain why in assemblies repulsive 
' spheres the boson peak shifts up to zero frequency as the pressure is lowered toward the jamming 

threshold, and derive the corresponding exponent. I show how these ideas capture the main low- 
frequency features of the vibrational spectrum of amorphous silica. These results extend arguments 
of Phillips^ on the presence of floppy modes in under-constrained covalent networks to glasses where 
r-{ ' the covalent network is rigid, or when interactions are purely radial. 
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Elasticity and transport in crystals are fairly well understood subjects, due to the simplicity of their underlying 
lattice. By contrast, in amorphous solids disorder has strong effects at intermediate length scales. As a consequence 
energy transport, force propagation and low-energy excitations are not yet satisfyingly understood. Another intriguing 
property is the ubiquitous presence in their vibrational spectrum of an excess of low-frequency modes with respect to 
the Debye prediction, the Boson Peak. Many observations support that this Peak is a key property of glasses: near 
the peak frequency transport is strongly affected^, the amplitude of the peak correlates well with the dynamics of the 
glass transition 6 and with the spatial extension of the force chains^. 

The cause and the nature of the peak are debated questions. Various approaches have focused on the role of 
disorder— i 10 ' 1:L i 12 ' 13 . In the present work we will follow a more geometric and perhaps more intuitive road, based on 
\ the concept of rigidity. Following ideas of Maxwell^, Phillips^ realized that as the composition of covalent glasses 
<— | ■ is changed to increase the mean valance of the atoms, the backbone of covalent bonds could undergo a rigidity 
transition^. Under-constrained covalent networks are floppy, and therefore present nearly zero-frequency modes 
O i (which get a small but non-zero frequency due to the presence of van der Waals interactions). Such floppy modes are 
indeed observed in the spectra of those weakly-connected solid s 16 i 17 , causing a boson peak whose amplitude tends to 
decrease as atoms of high valence are added to the composition of the glass, as predicted. These arguments explain the 
nature of the Boson Peak in under-coordinated covalent networks. It nevertheless leaves two fundamental questions 
unanswered: (i) why do glasses with a rigid, sufficiently connected backbone can still present a large Boson Peak, as 
is the case for example for silica? (ii) Can this description be applied to non-covalent glasses, known to present a 
JO I Boson Peak as well? In what follows I shall review recent results^ addressing these points. 

. ■ One inherent difficulty in the study of amorphous solids is that the length scales where the effects of disorder 
become strong is moderate. This lack of large parameter makes it hard to test in a stringent manner and distinguish 
clearly the consequences of different theories. It was proposed by Alexander— that the situation may be different in 
assemblies of purely repulsive particles, such as emulsions or elastic grains, as the pressure vanishes toward zero. This 
occurs at a packing fraction (j) c where repulsive particles are just in contact, called the "jamming transition", which 
corresponds to the so-called "random close packing" for mono-disperse spheres. He proposed that <f> c is a critical 
point. This idea was latter substantiated by the findings that the elastic modul i 19 ' 20 , the vibrational spectrum^ and 
force propagation 7 display scaling behaviors near the jamming threshold. As we shall see, although the system is 
amorphous and isotropic, it cannot be described as a continuous elastic body on any length scale. Because the strong 
effects of disorder occur already at large length scales near this critical point, this model system is a lens allowing to 
probe the properties of amorphous solids, and the effects of disorder. 

We shall focus on the vibrational spectrum, and its relation with the coordination of the packing. Numerical 
experiments 19 showed that at the jamming threshold, the spectrum of vibrational modes D(u>) does not present a 
Debye behavior D(u>) ~ lo 2 on any frequency range, rather one finds that D(uj) ~ u . There are no frequency-range 
were plane waves can be found. As the system is compressed, this "plateau" erodes below some frequency scale to*, 
see Fig([T]). In what follows we shall explain these behaviors and relate u>* to the coordination z of the packing. We 
shall then show how these concepts apply to silica. 
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FIG. 1: The density of vibrational states, D(lu), vs angular frequency, lu, for the simulation of Ref.— . 1024 spheres interacting 
with repulsive harmonic potentials were compressed in a periodic cubic box to packing fraction <j>, slightly above the jamming 
threshold <f> c . Then the energy for arbitrary small displacements was calculated and the dynamical matrix inferred. The curve 
labeled a is at a relative packing fraction cf> — (f> c = 0.1. Proceeding to the left the curves have relative volume fractions 10 -2 , 
10" 3 , 10~ 4 , 10" 8 , respectively. 

II. VIBRATIONAL SPECTRUM NEAR THE JAMMING THRESHOLD 

A. Energy expansion and soft modes 

Following 1 - we consider N soft spheres packed into a spatially periodic cubic cell of side L at volume fraction <fi. 
In what follows we consider repulsive, finite-range "soft spheres" . For inter-particle distance r < a, the particles have 
non-zero mutual energy and are said to be in contact. They interact with the following potential: 

V(r) = ~(l-~Y (1) 
a \ a) 

where tr is the particle diameter and e a characteristic energy. For r > a the potential vanishes and particles do not 
interact. Henceforth we express all distances in units of a, all energies in units of e, and all masses in units of the 
particle mass, m. In the following, we consider the harmonic case a = 2, but our arguments can be applied to other 
potential (for example Hertzian contacts where a = 5/2) 2 . In the harmonic case the energy expansion follows 2 -: 

^=te(^-i) [(4 :f )112 ] 

+^[(^-^).n y ] 2 +0(^ 3 ) (2) 

where the sum is over all iV c contacts (ij), r^J is the equilibrium distance between particles i and j, is the unit 
vector along the direction ij, and (SRj — SRi) 1 ^ indicates the projection of SRj — SRi on the plane orthogonal to n^. 
Note that Eq. ^ can be written at first order in matrix form, by defining the set of displacements 8Ri...5Rn as 
a (i/V-component vector |<5R). Then Eq. ([2]) can be written as SE = (SH\Ai\SH) . The corresponding matrix A4 is 
known as the dynamical matrix 21 -. The dN eigenvectors of the dynamical matrix are the normal vibrational modes 
of the particle system, and its eigenvalues are the squared angular frequencies of these modes. 

The first term in Eq. |2]) is proportional to the contact forces. Near the jamming transition r^J — > 1 so that this 
term becomes arbitrarily small. It has, nevertheless, interesting consequence on the mechanical stability of the system 
for <f> > <p c , and this term is important to understand how the coordination z varies with <f>, as discussed in^. In 
what follows we shall neglect it. This approximation corresponds to a real physical system where the soft spheres are 
replaced by point particles interacting with relaxed springs. We now have: 



5E = -Y}{5R 3 -5Ri).rl l3 f 



(3) 



FIG. 2: (color online) Illustration of the boundary contact removal process described in the text. Eighteen particles are confined 
in a square box of side L periodically continued horizontally and vertically. An isostatic packing requires 33 contacts in this 
two-dimensional system. An arbitrarily drawn vertical line divides the system. A contact is removed wherever the line separates 
the contact from the center of a particle. Twenty-eight small triangles mark the intact contacts; removed contacts are shown 
by the five white circles. 

If the system has too few contacts, there is a set of displacements modes of vanishing restoring force and thus vanishing 
vibrational frequency. These are the soft modes, or floppy modes. For these soft modes the energy 5E = of Eq. ([3|) 
must vanish; therefore they must satisfy the N c constraint equations: 

(SRi — 8Rj).fiij = for all N c contacts (ij) (4) 

This linear equation defines the vector space of displacement fields that conserve the distances at first order between 
particles in contact. The particles can yield without restoring force if their displacements lie in this vector space. 
Eq. (H]) is purely geometrical and does not depend on the interaction potential. Each equation restricts the dN- 
dimensional space of \5H) by one dimension. Except when specific packing symmetries are present ( which is not the 
case here), these dimensions are independent, so that the number of independent soft modes is dN — N c 40 . 

B. Isostaticity 

There are no internal soft modes in the rigid structure of a solid. This is true for a system of repulsive spheres, 
as soon as it jams 4 ^. Therefore jammed states must satisfy N c > dN, which is the Maxwell criterion for rigidity. At 
the jamming transition, this inequality becomes an equality, as was verified ir>22. Such a system is called isostatic. 
The coordination number z is then z c = 2N C /N — > 2d. Physically the reason why z — z c rather than z > z c is the 
following: if z > z c , the system is over-constrained and contacts must be straine d 23 ' 24 ! 25 . This cannot be so at the 
jamming transition where the pressure and all contact forces vanish. Thus, the jamming transition is very different 
from rigidity percolation models where springs are deposited randomly on a lattice. In these models, which are in a 
sense at infinite temperature, there are both over-constrained and floppy regions at the rigidity threshold, and the 
percolating rigid cluster is a fractal object with dimension smaller than d. In sphere packing the rigid system is a 
d-dimensional object. In silica, we shall see that the spatial fluctuations of coordination are small, and we will argue 
that the vibrations of this system resembles those of sphere packing. In chalcogenide glasses where atoms of different 
valence are used, the properties of the covalent network near the rigidity threshold are still debated, but it has been 
proposed that isostaticity may also characterize the rigidity threshold 2 ^. 

An isostatic system is marginally stable: if q contacts are cut, a space of soft modes of dimension q appears. For our 
coming argument we need to discuss the extended character of these modes. In general when only one contact (ij) is 



FIG. 3: One soft mode in two dimensions for N « 1000 particles. Each particle is represented by a dot. The relative 
displacement of the soft mode is represented by a line segment extending from the dot. The mode was created from a 
previously prepared isostatic configuration, periodic in both directions, following^. 20 contacts along the vertical edges were 
then removed and the soft modes determined. The mode pictured is an arbitrary linear combination of these modes. 

cut in an isostatic system, the corresponding soft mode is not localized near (ij). This arises from the non-locality 
of the isostatic condition that gives rise to the soft modes, see^ for a theoretical derivation. When many contacts 
are severed, the extended character of the soft modes that appear depends on the geometry of the region being cut. 
If this region is compact many of the soft modes are localized. For example cutting all the contacts inside a sphere 
totally disconnects each particle within the sphere. Most of the soft modes are then the individual translations of 
these particles and are not extended throughout the system. In what follows we will be particularly interested in the 
case where the region of the cut is a hyperplane as illustrated in Fig. ^j. In this situation occasionally particles in the 
vicinity of the hyperplane can be left with less than d contacts, so that trivial localized soft modes can also appear. 
However we expect that there is a non-vanishing fraction q' of the total soft modes that are not localized near the 
hyperplane, but rather extend over the entire system, like the mode shown in Fig. ([3]). This assumption is checked 
numerically irA We shall define extended modes more precisely in the next section. 

C. D(lj) of an isostatic configuration 

We aim to show first that the density of states of an isostatic system does not vanish at zero frequency. D(u>) is the 
total number of modes per unit volume per unit frequency range. Therefore we have to show that there are at least 
on the order of ujL d normal modes with frequencies smaller than uj for any small ui in a system of linear size L. As we 
justify later, if proven in a system of size L for lu ~ lul ~ 1/L, this property can be extended to a larger range of u) 
independent of L. Therefore it is sufficient to show that they are of the order of L d ~ x normal modes with frequency 
of the order of 1/L, instead of the order of one such mode in a continuous solid. Our procedure for identifying the 
lowest frequency modes resembles that used for an ordinary solid. An isolated block of solid has three soft modes 
that are simply translations along the three co-ordinate axes. If the block is enclosed in a rigid container, translation 
is no longer a soft mode. However, one may recover the lowest-frequency, fundamental modes by making a smooth, 
sinusoidal distortion of the original soft modes. We follow an analogous procedure to find the fundamental modes of 
our isostatic system. First we identify the soft modes associated with the boundary constraints by removing these 
constraints. Next we find a smooth, sinusoidal distortion of these modes that allows us to restore these constraints. 

For concreteness we consider the three-dimensional cubical TV-particle system S of Ref. 19 with periodic boundary 
conditions at the jamming threshold. We label the axes of the cube by x, y, z. S is isostatic, so that the removal of 
n contacts allows exactly n displacement modes with no restoring force. Consider for example the system S' built 
from S by removing the q ~ L 2 contacts crossing an arbitrary plane orthogonal to (ox); by convention at x — 0, see 
Fig. S 1 , which has a free boundary condition instead of periodic ones along (ox), contains a space of soft modes of 
dimension q^, instead of one such mode — the translation of the whole system — in a normal solid. As stated above, 
we suppose that a subspace of dimension q' ~ L 2 of these soft modes contains only extended modes. We define the 
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extension of a mode relative to the cut hyperplane in terms of the amplitudes of the mode at distance x from this 
hyperplane. Specifically the extension e of a normalized mode |<5R) is defined by sin 2 (^i^)(i|<5R) 2 = e, where the 
notation (i|JR) indicates the displacement of the particle i of the mode considered. For example, a uniform mode 
with (i|5R) constant for all sites has e = 1/2 independent of L. On the other hand, if (i|<5R) = except for a site i 
adjacent to the cut hyperplane, the Xi/L ~ L -1 and e ~ L~ 2 . We define the subspace of extended modes by setting 
a fixed threshold of extension eo of order 1 and thus including only soft modes (3 for which ep > cq. As we discussed 
in the last section, we expect that a fixed fraction of the soft modes remain extended as the system becomes large. 
Thus if q' is the dimension of the extended modes vector space, we shall suppose that q' /q remains finite as L — > oo. 
The appendix presents our numerical evidence for this behavior. 

We now use the vector space of dimension q 1 ~ L 2 of extended soft modes of S 1 to build q' orthonormal trial 
modes of S of frequency of the order 1/L. Let us define |<5R,g) to be a normalized basis of this space, 1 < /3 < q'. 
These modes are not soft in the jammed system S since they deform the previous q contacts located near x = 0. 
Nevertheless a set of trial modes, \SR*p), can still be formed by altering the soft modes so that they do not have an 
appreciable amplitude at the boundary where the contacts were severed. We seek to alter the soft mode to minimize 
the distortion at the severed contacts while minimizing the distortion elsewhere. Accordingly, for each soft mode (3 
we define the corresponding trial-mode displacement (i|<5R*) to be: 

(^R^C^sm^X^) (5) 

where the constants Cp are introduced to normalize the modes. Cp depends of the spatial distribution of the mode 
p. If for example, a highly localized mode has (i\SH) = except for a site i adjacent to the cut plane, Cp grows 
without bound as L — > oo. In the case of extended modes C^ 2 = Yluj) sm2 ( £ |r)(j : \SRp ) 2 = e p > e o, and therefore 

— 1/2 

Cp is bounded above by e . The sine factor suppresses the problematic gaps and overlaps at the q contacts near 
x = and x — L. The unit basis I^R^) can always be chosen such that the |(>R^) are orthogonal, simply because 
the modulation by a sine that relates the two sets is an invertible linear mapping in the subspace of extended modes. 
Furthermore one readily verifies that the energy of each |<5Rp is small, and that the sine modulation generates an 
energy of order 1/L 2 as expected. Indeed we have from Eq. ^ : 

SE = C 2 ^[(sin^Xzl^) - S in(^)(j|^)) • ft^f (6) 

m 

Using Eq. ((4]), and expanding the sine, one obtains: 

SE « C 2 ]T cos 2 (^)^(n y - • e x ) 2 ((j\SRp) ■ n l0 f (7) 

m 

<e \n/L) 2 J2(j\S^) 2 (8) 

m 

where e x is the unit vector along (ox), and where we used | cos | < 1. The sum on the contacts can be written as a 
sum on all the particles since only one index is present in each term. Using the normalization of the mode (3 and the 
fact that the coordination number of a sphere is bounded by a constant z max (z max — 12 for 3 dimensional spheres^), 
one obtains: 

SE < e^^/Lfz^ = uj\ (9) 

We have found on the order of L 2 trial orthonormal modes of frequency bounded by wl ~ 1/L. This leads to 
D(u>l) ~ L 2 /lul ~ L°, i.e. the average density of states is bounded below by a constant below frequencies of the 
order u>l. This scaling argument is not entirely rigorous as the orthonormal modes we use are not the normal modes. 
For completeness, we show how a variational procedure leads to a rigorous proof. M. is a positive symmetric matrix. 
Therefore if a normalized mode has an energy SE, we know that the lowest eigenmode has a frequency ujq < \J SE. 
Such argument can be extended to a set of modes. If m a is the a'th lowest eigenvalue of M. and if e a is an orthonormal 
basis such that (e a \M\e a ) = n a then the variational bound of A. Horn [Am. J. Math 76 620 (1954)] shows that 
Yl m a < Yl n a- Since 1 n q ^ S^o, and since YX m a > J2 q q /2 m a ^ (<l/' 2 )m q /2, we have qm p > {q/2)n q / 2 as 
claimed. Thus, if there are m orthonormal trial modes with energy SE < us 2 , then there are at least m/2 eigenmodes 
with frequency smaller than y2w ( . Therefore finding of the order of L d ~ x trial orthonormal modes with energy of 
order 1/L 2 indeed leads to a non- vanishing density of normal modes. 
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FIG. 4: (color online) Scaling of u* with the excess coordination number, Sz in the system with relaxed springs. The line has 
a slope one. 

In what follows, the trial modes introduced in Eq. ([5]), which are the soft modes modulated by a sine wave, shall 
be called "anomalous modes" to distinguish them from plane waves. 

To conclude, one may ask if this variational argument can be improved, for example by considering geometries of 
broken contacts different from the hyperplane surfaces we have considered so far. When contacts are cut to create 
a vector space of extended soft modes, the soft modes must be modulated with a function that vanishes where the 
contacts are broken in order to obtain trial modes of low energy. On the one hand, cutting many contacts increases 
the number of trial modes. On the other hand, if too many contacts are broken, the modulating function must have 
many "nodes" where it vanishes. Consequently this function displays larger gradients and the energies of the trial 
modes increase. Cutting a surface (or many surfaces, as we shall discuss below) is the best compromise between these 
two opposit effects. Thus our argument gives a natural limit to the number of low-frequency states to be expected. 

We may extend this argument to show that the bound on the average density of states extends to higher frequencies. 
If the cubic simulation box were now divided into m 3 sub-cubes of size L /m, each sub-cube must have a density of 
states equal to the same D(u>) as was derived above, but extending to frequencies on order of mujL. These subsystem 
modes must be present in the full system as well, therefore the bound on D{uj) extends to [0,towl]. We thus prove 
that the same bound on the average density of states holds down to sizes of the order of a few particles, corresponding 
to frequencies independent of L. We note that in d dimensions this argument may be repeated to yield a total number 
of modes, L d_1 , below a frequency ujl ~ thus yielding a limiting nonzero density of states in any dimension. We 
note that the trial modes of energy SE ~ Z _1 that we introduce by cutting the full system into subsystems of size I 
are, by construction, localized to a distance scale I. Nevertheless we expect that these trial modes will hybridize with 
the trial modes of other, neighboring, subsystems; the corresponding normal modes will therefore not to be localized 
to such short length scales. 



When the system is compressed and moves away from the jamming transition, the simulations show that the extra- 
coordination number Sz = z — z c increases. This causes AN C = NSz/2 ~ L d Sz extra constraints to appear in Eq. ([3]). 
Cutting the boundaries of the system, as we did above, relaxes q ~ constraints. For a large system, L d Sz > L' 1 ^ 1 
and thus q < AN C . Thus Eq. ([3]) is still over-constrained and there will be no soft modes in the system. However, as 
the systems become smaller, the excess number of constgraints, AN C , diminishes; for L smaller than some I* ~ 5z , 
q becomes larger than A7V c the system is again under-constrained as was already noticed in^ 3 -. This allows one to 
build low-frequency modes in subsystems smaller than I*. These modes appear above a cut-off frequency uj* ~ Z* -1 ; 
they are the "anomalous modes" that contribute to the flat plateau in D(oj) above u>* . In other words, anomalous 
modes with characteristic length smaller than I* are not affected very much by the extra contacts, and the density of 
states is unperturbed above a frequency lu* ~ Sz. This scaling is checked numerically in Fig. |4j This prediction is in 
very good agreement with the data up to Sz ~ 2. 

At frequencies lower than lo* we expect the system to behave as a disordered, but not poorly-connected, elastic 
medium. These vibrational modes are similar to the plane waves of a continuous elastic body. We refer to these 
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modes as "acoustic modes" . The behavior of such systems near the jamming threshold thus depend on the frequency 
lo at which they are observed. For ui > u* the system behaves as an isostatic system, and for u> < u>* it behaves as a 
continuous elastic medium. Equivalently I* is the distance below which a continuous elastic description is not a good 
approximation. This was confirmed numerically where it was found that I* characterize the stress heterogeneities 
after a point force is imposed. 



III. THE CASE OF AMORPHOUS SILICA 



Silica is perhaps the most common glass. It is also known to have one of the strongest excess of low-frequency 
modes, or boson peak, see^I for a review of empirical results and models. In this paragraph we propose an explanation 
for its density of states at low-frequency, and for the nature of the excess-modes. In this glass (or more generally 
aluminosilicates) the forces within the tetrahedra SiO^ are much stronger than the forces that act between them 3 -^: 
it is easier to rotate two linked tetrahedra, that is to bend the Si-O-Si angle, than to distort one tetrahedron: the 
bending energy of Si-o-Si is roughly 10 times smaller than the stretching of the contact Si-o22. This suggests to 
model such glass as an assembly of linked tetrahedra loosely connected at corners: this is the "rigid unit modes" 
model 30 . In this model the tetrahedra are characterized by a unique parameter, a stiffness fc 44 . Recently this model 
was used to study the vibrations of silica^. The authors first generate realistic configurations of S1O2 at different 
pressures using molecular dynamics simulations. At low pressure, they obtain a perfect tetrahedral network. When 
the pressure becomes large, the coordination of the system increases with the formation of 5-fold defects. Once these 
microscopic configurations are obtained, the rigid unit model is used and the system is modeled as an assembly of 
elastic tetrahedra, see Fig.©. Then, the density of states of such network is computed. The results are shown in 
Fig.©. One can note the obvious similarity with the density of states near jamming of Fig.}]}. We argue that the 
cause is identical, and that the excess-modes correspond to the anomalous modes made from the soft modes, rather 
than to one-dimensional modes as proposed in2&. Indeed, a tetrahedral network is isostatic, see e.gi^l. The counting 
of degrees of freedom goes as follows: on the one hand each tetrahedron has 6 degrees of freedom (3 rotations and 
3 translations). On the other hand, the 4 corners of a tetrahedron bring each 3 constraints shared by 2 tetrahedra, 
leading to 6 constraints per tetrahedron. Thus the system is isostatic. When the pressure increases the coordination 
increases too, leading to the erosion of the plateau in the density of states discussed earlier. 

These predictions fail to describe the spectrum of silica at low frequencies, where the weak interactions, in particular 
the bending of the Si-O-Si bond and the Van der Waals interactions cannot be neglected. The main effect of these 
interactions is to shift the spectrum of anomalous modes by some frequency scale £■ To estimate Q we use the stiffness 
of the Si-O-Si bending interaction obtained ab initio^, and the molecular mass to form a frequency. One finds 
£ « lAThz. We therefore expect the spectrum of silica to display a plateau above a frequency of order 1 Thz. This 
is indeed what is observed in simulations: silica glass present a well-defined plateau above the boson peak frequency, 
as appears in the numerical results shown in Fig. ([7]). 

Our argument does not involve disorder. Thus it must also apply for the crystals of the same composition and similar 
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FIG. 6: Density of rigid unit modes for silicate at different pressure. Trachenko et al.— . 



densities such as the a and /3-cristobalite, since these crystalline structures are formed, as silica, by SiC>4 tetrahedra 
connected at the corners. /3-cristobalite has the structure of the diamond, in which the tetrahedra correspond to the 
4 carbons bonded to a central carbon, whereas a-cristobalite has a tetragonal structure. Empirically a boson peak is 
observed in all these materials 53 . Numerically, a plateau indeed appears in D(to) at roughly the same frequency in 
the cristobalite a and f3^1 and in the glass, as shown in Fig. ([7]). In the crystalline case the plateau corresponds to a 
Van Hove singularity and to the accumulation of optical bands. More generally, there are other crystals showing an 
excess density of states at frequencies of the order of the boson peak frequency of the corresponding glas o 31 ' 34 i 35 ' 36 . 
According to the present argument this should be the case as long as the crystal and the amorphous structure share 
a similar connectivity (which is not the case, for example, for sphere packing). 

Although disorder is not relevant to compute the density of states, it affects the nature of the vibrational modes. 
Transport is very different in silica glass and in cristobalite. Thus the peculiarity of the amorphous state lies in the 
nature of the excess-modes, not in the density of states 3 - 1 -. It is useful to note the parallel between cristobalite and 
silica glass on the one hand and cubic lattice and the jamming threshold of elastic spheres on the other hand. In both 
cases the amorphous solid and the crystal have a similar density of states, but the anomalous modes in the amorphous 
phase are not plane waves. Disorder strongly affects the anomalous modes, and makes them very heterogeneous as 
shown in Fig $5§. 

In chalcogenide glasses, assuming weak fluctuations of coordination (which may be a naive assumption in that case) 
our argument leads to a Boson Peak frequency increasing for over-constrained networks as co* ~ (£ 2 + a(x — a^) 2 ) 1 / 2 , 
where x c is the critical composition above which floppy modes vanish, a is a numerical constant and £ is the finite 
frequency of the floppy modes for x < x c , induced by van der Waals interactions. When the coordination becomes 
too large the Boson Peak vanishes, as occurs for example for amorphous silicon. 



IV. CONCLUSION 



There exists several analytical theories of the Boson Peak based on specific models of the disorder. Some of the 
authors have investigated models in which the disorder enters via local defects 3 -!; intrinsic disorder in the harmonic 
force constants, see e.gj&i^, in which the boson peak frequency marks a crossover from acoustic-like sound excitations 
to a disorder-dominated regime, or infinite temperature systems where particle positions are rando m 11 ! 12 ' 13 . These 
models are able to generate a Boson Peak, to explain the presence of a plateau in the thermal conductivity of glasses, 
see e.g£, and can describe sound attenuation at those frequencies^. With respect to those approaches, the rigidity, 
or coordination -based description has the advantages of: (i) relating the boson peak to a microscopic observable, 
the coordination number, which can be inferred from the composition of the glass for covalent networks, and more 
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FIG. 7: Density of states of silica glass (at temperatures of 10 and 300 Kelvins), a-cristobalite and /3-cristobalite. This figure 
is taken from the simulations o^. 



generally from the knowledge of the microscopic structure (ii) explaining why crystalline structures can also present 
a peak and allows to compute its position (iii) introducing a length scale I* characterizing the modes responsible for 
the peak. This length scale was shown to describe the extension of the force chains^. These results can be extended 
to colloidal glasses^, and systems with rapidly decaying long-range potentials, such as a Lennard- Jone o 3 i 39 . 
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Of these, d(d + l)/2 modes are dictated by the translational and rotational invariance of the energy function <5_E[<5R]. Apart 
from these, there are dN — N c — d(d + 1) /2 independent internal soft modes. In what follows we shall neglect the term 
d(d+ l)/2, as it effect vanishes for large N. 

The rattlers, particles without contacts, are removed in this argument, as they will lead to trivial zero translational modes. 
The balance of force can be satisfied in S' by imposing external forces on the free boundary. This adds a linear term in the 
energy expansion that does not affect the normal modes. 

In a polydisperse system z max could a priori be larger. Nevertheless Eq. (|8]) is a sum on every contact where the displacement 
of only one of the two particles appears in each term of the sum. The corresponding particle can be chosen arbitrarily. It 
is convenient to choose the smallest particle of each contact. Thus when this sum on every contact is written as a sum on 
every particle to obtain Eq. ©, the constant z max still corresponds to the monodisperse case, as a particle cannot have 
more contacts than that with particles larger than itself. 

In fact the rigidity of a tetrahedron induced by the covalent bonds should be characterized by 3 parameters corresponding 
to different deformations of the tetrahedron. If these parameters are of similar magnitude, as one expects for example for 
silica, this does not change qualitatively the results discussed here. 
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